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Abstract 



Q I As a model of coupled nano-electromechanical resonantors we study two nonlinear driven os- 

' cillators with an arbitrary coupling strength between them. Analytical expressions are derived 

for the oscillation amplitudes as a function of the driving frequency and for the energy transfer 
rate between the two oscillators. The nonlinear restoring forces induce the expected nonlinear reso- 
nance structures in the amplitude-frequency characteristics with asymmetric resonance peaks. The 
' corresponding multistable behavior is shown to be an efficient tool to control the energy transfer 

\ that the nonlinear response can be exploited to design precise sensors for mass or force detection 



arising from the sensitive response to small changes in the driving frequency. Our results imply 



experiments based on nano-electromechanical resonators. 
PACS numbers: 



I. INTRODUCTION 



Currently experimental efforts are devoted to the fabrication of nanoscale resonators with 
a precise control of their behavior. Such nanoscale resonators are ideal prototype systems for 
testing fundamental physical concepts, such as entanglement and quantum correlations 
By now several types of resonators were successfully considered, like optical two level atoms 

n n n 

in quantum cavities j2|, artificial Josephson junction qubits [3|, atoms seized in ion traps 
or nano-optomechanical devises l5|. There is a new trend towards nano-electromechanical 

n n 

resonators. They are widely studied from both the quantum-mechanical |6| and from a 



classical point of view |7H14|. These devices are approximately 200 nm in size and consist 
of three layers of gallium arsenide (GaAs): an n-doped layer of width 100 nm is stacked 
within an insulating layer of 50 nm and a p-doped layer of 50 nm [l^. The resonators 
can be controlled by electric fields via the piezoelectric effect, which fix their mechanical 



strain 



13| . Along with a single resonator, one can consider coupled resonators driven by an 



additional external field. Coupled resonators show different dynamical regimes dependent on 
the interplay of their coupling strength and the driving. For the case of moderate coupling 
between two resonators this problem was already addressed in a recent paper by Karabalin 
et al. [l^. They showed that the linear and weakly nonlinear response of one oscillator 
can be modified by driving the second oscillator. A complicated frequency- sweep response 
curve was obtained numerically when both oscillators are driven into the strongly nonlinear 
regime. 

In this paper we study two nonlinear oscillators allowing for an arbitrary coupling strength 
between them with a possibility of driving both with the same frequency but different am- 
plitudes. The coupling strength between the oscillators is quantified in terms of the connec- 
tivity parameter defined below in Sect. II. We derive general analytical expressions for the 
amplitude-frequency characteristics valid for arbitrary (weak as well as strong) connectivity. 
We analyze the redistribution of energy between the two resonators injected into the system 
via the external driving fields. We quantify stable and unstable dynamical regimes, with 
special focus on the nonlinear response of the system. Our predictions point to possible 
new applications of nanoscale resonators exploiting their sensitivity in response to external 
: ields and perturbations. In particular they may be used as sensors for tiny forces or masses 



15 



16| which lead to a shift in their resonance frequencies to be identified in the sensitive 



2 



nonlinear response regime. 

Our paper is organized as follows: in the next section, we discuss our fundamental model 
of two coupled driven oscillators. In Sect. Ill we study the mode frequency shifts and 
relaxation effects for a non-resonant driving, while in Sect. IV we address the resonant 
case with a special focus on the nonlinear shifts of the mode frequencies. In the following 
sections we investigate the frequency response function and the key problem of the energy 
redistribution between the oscillators, before concluding in Sect. VIII. 



II. MODEL 



We consider two nanomechanical oscillators described by the coordinates xi^2 in the frame- 
work of the model outlined in [l^. The corresponding dynamical equations can be written 
down in the following form: 

Xi + culxi + D{xi — X2) = eM, 

X2 + 00^X2 + D{x2 - xi) = eN. (1) 
eM = eM(xi, xi, t) = —271X1 — xi^l + Fi cos fit, 
eN = eN{x2, ±2, t) = -272X2 - X2xl + F2 cos fit, (l') 

£ < 1. 

Here cui and U2 are frequencies of the individual resonators, 71 and 72 are the dissipation 
coefficients, Xi ^-^d X2 are nonlinearity parameters, Fi and F2 are amplitudes of the external 
harmonic forces applied to the resonators, fi is the frequency of these forces and D is the 
coefficient of the linear coupling between the resonators. As usual, we assume the right hand 



side of Eqs. (1) to be small perturbations, see also 

We summarize the canonical solution of the unperturbed coupled system first. Its dy- 
namics follows the simple equations 

xi + u^xi + D{xi - X2) = 0, 

X2 + (^2^2 + -D(X2 - Xi) = 0. (2) 

The transition from coupled oscillations to the mode oscillations can be done via the following 
transformation 20| . 

xi = qi + q2, 



X2 = -K-\i + Kq2, (3) 

where 

K = --(l + v^IT^), K-^ = -{l-VlT^), KK-^ = 1. (4) 
a a 

2D , . 

^ = VT^-TTr (5) 

We call the parameter a describing the coupling strength between the oscillators connectivity. 
From now on we assume that UJ2 > Wi- 

The mode oscillations have the frequencies 



z/?o =w^Ta;Vl + a2, (6) 



'1,2 

where 

~2 _ + ^2 
^1,2 = ^1,2 + ^) 

(7) 

0)1^2 are partial frequencies. We would like to stress that the value of the connectivity a 
depends not only on the linear coupling term D, but on the proximity of the free oscillation 
frequencies Ui and 002- In the limit of a weak connectivity (a ^ 1) the frequencies z/i_2 tend 
to the partial frequencies 0)1,2, while in the limit of strong connectivity (a ^ 1), 

~ ul, vl ~ (S^, (8) 

where 



Cu I — , CU I — -v . 



2 ' ^ 2 

Obviously, in the limit a ^ 1 the mode frequency separation attains the maximal possible 
value which is equal to 2D. 

In the case of a finite driving Fi 2 7^ but in the linear (xi,2 = 0) dissipationless regime 
(71,2 = 0) the particular solutions of the dynamical equations (1) are given by 

Xi = Ai cos Qt, X2 = A2 cos Qt, 

F,,2i^l - + F2,lD .QN 

Ai,2 J, , (9) 
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here 

1 _ 1 

d2 ~ (l]2_;,2)(^2_j,2) 



I \ 2 \ ^^Q^ 



and 2 are the amphtudes of the induced resonator oscillations. They increase resonantly 
when the frequency of the external driving tends closer to either of vi^2- 

With the above solutions at hand the influence of the damping term as well as of the 
nonlinear corrections can be taken into account with the help of a standard substitution in 
the resonant denominator [2^: 

1 1 1 

-, :rT : z = ^ -)■ 



1,2 



9? 2z/i^2('^l,2 - fi) 2z/i^2('^l,2 + ^1,2 + ^ri^2 - fi) 



; (11) 

2i^l,2JK2 + 51,2-^^)2 + ^2 



where Fi 2 are the mode relaxation rates, 81^2 are the nonlinear corrections to the mode 
frequencies that depend on the oscillation amplitudes A1.2. We would like to point out that 
the substitutions of Eq. (fTTj) are correct only if the nonlinearity is not too strong and the 
decay rate is not too high (i^i 2 ^ ^1,2; ri,2)- Explicit expressions for 81^2 and Fi 2, in terms 
of the system parameters will be given in the next section. 

III. NONLINEAR SHIFT OF THE MODE FREQUENCIES AND CONSE- 
QUENCES OF THE RELAXATION TERMS. THE NON-RESONANT CASE 

We turn back to the perturbed system of Eq. (1) making the preliminary assumption 
that the resonance condition does not hold at the mode frequencies 1^1,2 7^ ^- It is clear that 
in this particular case, the role of the external force is negligible. We concentrate at first 
on the influence of the relaxation and of the nonlinearity on the oscillation of the coupled 
resonators. 

To study the equations (1), we use a modified method of a slowly varying amplitudes, 
[itI . Taking into consideration the transformation (3), we reduce the unperturbed system 
to the mode oscillations and write down the solution in the form: 

xi{t) = Ai{t) sm[uit + ai{t)] + A2{t) sm[u2t + a2{t)], 
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X2it) = -K-^Ai{t) sin[z/it + ai{t)] + KA2{t) sin[z/2t + aaW], (12) 

where Ai^2(t), ai,2(^) are slowly varying amplitudes and phases, respectively. The variables 
Ai 2{t) are the first-order infinitesimal variables and therefore the terms proportional to the 
second-order derivatives Ai 2(t) can be omitted upon an insertion of the above relations into 
Eq. (1). Following the standard procedures, after straightforward but laborious calculations, 
we find for the slowly varying amplitudes and the phases 

:{KP2+Q2), 



dt Avi Vl + 
dA2 1 cr 



dt 4i/2 yjl + a^ 

^ dt 41/1 v/fT^^ ' 

= ^^T^iKP, + Q,), (13) 
dt AV2 Vl + 



where 



■j^ /»27r /'27r nliv /'27r /»27r /'27r 

Pl = ^J^ M COS idid7]dC Qi = ^J^ Ncos^d^d7]dC, 

2 p27v p27r p27r ^ p27r p27v p27r 

P2 = — ^/ / / M COS r]dCdr]dC, Q2 = -r^ / / N cosrjd^dridC, 



4vr3 7o J, J, ' ^ ' ^' 47r3 



^0 ^0 



-1 p27r p27v p27r -1 p27T p27r p27r 

P3 = -r^ / / MsinCdCdridC Qs = — / / N sin ^d^drjdC 



47r3 7o 7o 7o - - ' - 47r3 



^0 JO 



-| p27v p27r p27r -1 p27r p27r p27r 

P4 = — / / / Msinrid^dridC, Q4 = — / / / N smrid^drjdC. (14) 
For M and A^, determined by Eq.(l'), after inserting Xi{t) and X2{t) from (fT2l) . we obtain 

r], = — 7i(y4ii/i cos.^ + A2i'2 cos 77 + AiQ sin ^) — ai{Ai sin^ + y42 sin r] + Ai cos ()^, 

N{^,VX) = -72(-^A^-VicosC + ^2^^2/2 COST/- ^2^ sin C)- 
-a2(-^iA'"^ sin^ + A2A'sinr7 + A2 cosC)^, 

^ = uit + ai; ?7 = z/2^ + a2; C = (15) 
Inserting (15) into (14), after simple integration, from (13) one infers 

dt 2^' dt 2 
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dai da2 



where 



1 

2 L 



dt 



71 1± 



dt 



52; 



(16) 



VTTa 



72 1 T 



3 1 r .2f- . 1 



1^1,2 



a' 



X2^2 It 



0"^ 



(17) 



are the relaxation rates and the nonhnear shifts of the mode frequencies. An interesting fact 
is that the shift of the mode frequencies 81^2 depends on the square of the amphtudes 2 
as a consequence of the nonhnearity. 

In the case of a weak connectivity (cr ^ 1), from Eq. (17) we deduce 

3 1 

ri,2 ^ 7l,2, ^1,2 = Xl,2^?,2, 



4wi : 



while in case of strong connectivity (cr ^ 1) 



ri = r2f«-(7i + 72), 



5i 



(18) 



-XX1AI + X2AI) 



Therefore, when cr ^ 1 the modes are damped with the equal rates. However, the 
nonlinear shifts of the mode frequencies are different (5i > 52). 



IV. FREQUENCY RESPONSE FUNCTION OF TWO COUPLED RESONATORS 



The amplitudes of the driven oscillations of the two coupled resonators are presented by 
the expressions (9) and (10). Reexpressed in terms of the connectivity cr, we can rewrite the 
amplitudes as 

i^l,2(c^2'l-^')+^2lwV 

^^'^ = '—W^ ^ ^ 

ujI ^ ' ' ' ' 



In the limit of a weak connectivity [o ^ 1) Ax^2 = 2-Fi.2/(i^i2 ~ ^^) whereas in the limit 



of a strong connectivity (cr ^ 1) we obtain 



DF2,i 



1 1 

+ 



1]2 



~,2 



(20) 



Nanomechanical resonators for intermediate values of the connectivity {a ^ 2) were 
studied numerically and experimentally in [l^. While the general analytical expression 
(19) is derived for arbitrary values of the connectivity and therefore in the special case of a 
moderate connectivity the solution recovers the amplitude-frequency characteristics obtained 
in [l^. Using the transformation (11), one can easily modify (19) in order to add corrections 
describing the damping and nonlinear terms. However, the expressions obtained in this way 
are rather involved (see Appendix A). That is why here we only present the asymptotic 
expressions corresponding to the strong connectivity limit (a ^ 1): 

where 61^2 and ri^2 are determined in (18). 

As follows from (21), in the case of a strong connectivity the force F2 acting on the second 
oscillator, "drives" the first one, and vice versa - Fi, acting on the first oscillator, "drives" 
the second one. 

The amplitude-frequency characteristic consists of two tilted peaks with different heights, 
see Fig. [H The first peak corresponds to the frequency u+ and is definitely more pronounced 
than the second peak, corresponding to the frequency a7+ (0/+ ^ u+). Furthermore, the 
first peak is more tilted due to the relation 61 > 62- The parts of the plot CD and IH, 
corresponding to unstable oscillations of the system, are dotted. During upward/downward 
frequency sweeps of Q one observes hysteretic behaviour around Aa;+ and Aa)_|_ along the 
loops -BCED and GIKH, respectively. In Q similar hysteretic loops in amplitude- 

frequency characteristics of coupled nonlinear oscillators were obtained numerically and were 
confirmed experimentally for intermediate values of the connectivity cr ^ 2. In the unstable 
region the system is extremely sensitive to the perturbations. This fact can be used for 
the switching of the oscillation amplitude. After reaching the point C, the amplitude of 
the oscillation decreases sharply to the value E. Therefore, a simple and efficient switching 
protocol can be realized by tuning of the external field frequency only. 

We would like to point out, that the domain of the amplitude frequency characteristics 
that should be utilized for switching belongs to the unstable area (see the frequency inter- 
vals B — C and G — / on Fig. 1). Therefore, the system can jump to the lower state before 
reaching the summit of the unstable domain (point C). If this happens the jump of the oscil- 
lation amplitude is smaller making difficult the experimental observation of the two different 
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FIG. 1: Amplitude- frequency characteristics for the system of two strongly coupled oscillators, 
plotted using Eq.(21) and following values of the parameters: Fi = F2, A = Ai^2, ^+ = 
1.07 • 10*^ Hz, = 1.2 • 10^ Hz, 61 = 50.4 • lO^M^ fj^^ ^ 45 . lO^M^ Hz, Ti = Ta = 
2.0 • 10^ Hz, £>Fi,2 = 81.6 • lO^^mHz^. 

transport regimes. To circumvent this problem the frequency of the driving field should be 
changed adiabatically. Following our approach here we seek the criteria of adiabaticity that 

ar problem arises for example 



2l|. The method of the slow 



may be useful for different realizations of the system. A simi 
when studying nonlinear resonant transport in cold atoms 
varying amplitudes implies that the amplitude change rate should be slower than the mode 
frequencies i^i^2- Therefore, the rate of the oscillation amplitude change caused by tuning 
the frequency of the driving field is limited by the following condition: 

dt dT < (22) 

The adiabaticity condition can then be simplified taking into account the explicit expressions 
for the amplitudes given in Eq. (21). It is not difficult to show that in the vicinity of unstable 
areas 

dn{t) ^ n{t) ■ ^ ^ 

Thus, for the adiabaticity criteria we finally obtain the following estimation 

^ < min(.,.)M. (24) 
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V. NONLINEAR SHIFT OF THE MODE FREQUENCIES AND THE INFLU- 
ENCE OF THE RELAXATION TERMS. THE RESONANT CASE 



Let us suppose that the harmonic force Fi cos f2t is tuned in resonance with one of the 
modes and that F2 = 0. For this problem we derive equations for the slowly varying 
amplitudes in a more straightforward way. Taking into consideration the resonance condition 
and the transformation (3), we can write down the solution of the equation set (1) in the 
following form: 

xi{t) = Ai{t) sinz/it + A2{t) cosi/it + B{t) sin(i/2t + ip{t)), 

X2{t) = -R-^ {Ai{t) sinuit + A2{t) cos uit + KB{t)sm{iy2t + ip{t)). (25) 

After application of the standard method outlined in the last Section, for equations of slowly 
varying amplitudes and phases we obtain 

^_ J ^^^A-ipW + gW), 



where 



dt 4z/i Vl + 



dt 4:1^1 y/l + a^ 

dt - Au2VTT^^ ' 



B dt Ap2 VTTa^ 



{KPr + QD, (26) 



Pi^ = ^j^ I M cos vd^dr], qW = ^^ ^ N cos vd^dT], 

= A r r M sin r^d^dv, Qt^ = r NsinT^d^dr], 

27r^ Jo Jo 2tx^ Jo Jq 



P^" = — / Mcos^d^dr], Q^;' = — / Ncosididiq, 
Jq Jo 27r^ Jo Jo 

Pt^ = 7^ r r M sin ^d^dv, = 7^ r r N sin ididrj. (27) 

27r^ Jo Jo 27r^ Jo Jo 

^ir) = M^'') (xi , xi , t) = -7i [i/i ( COS r/ - sin r/) + 5 cos - 
—Xi{Ai sin?7 + A2 cos?7 + B cos^)'^ + F cos fit, 
jyM _ ArM(a,2,i;2,t) = -72[-K-Vi(Ai cosr/ - sinr/) + z/a^cos^]- 
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- X2 [-K-^ {Ai sin r] + A2 cos rj) + KB sin (28) 

Upon insertion of Eqs.(28) into (27) and after an integration, one gets equations for the 
slowly varying amplitudes Ai,A2,B and ip. The explicit form of these equations is given in 
Appendix B. In the limit of a strong connectivity (cr ^ 1), the expressions simplify to 

As is evident, a resonant external force Fi cos fit, {Q = ui), for 7 = x = leads to the 
simplest form of instability (the secular instability), namely to the linear growth of the 
oscillation amplitude Ai = {F/4Q)t. 

We would like to point out that in the first three equations for the variables Ai, A2 and 
B, the right hand side of the set of Eqs. (29) does not depend on the fourth variable ip. 
Therefore, the set of Eqs. (29) can be solved self-consistently for the first three variables. 

In order to find the stationary values of the slowly varying amplitudes and in order to 
examine the stability of these values, we utilize the following transformation: 

Ai = pcos6', v42 = -psin6'. (30) 

In the more convenient polar coordinates p and 6 we obtain 

- = -7P + ^cos^, 

de _ F sine 

It - ~ 4fi^' 

where Uj^l = |^(p^ + 25). By setting the rhs of Eqs. (31) equal to zero, we obtain equations 
for the stationary values of amplitudes: 

i?o = 0, scos6'o = po5 ■ssin6'o = — rpQ , (32) 

where s = = f To determine po, we eliminate the variable Oq from the set of 

Eqs. (31) and obtain a cubic equation with respect to x = p^: 

-' + ^-^ = 0^ (33) 
11 



Eq. (33) is a reduced cubic equation. The number of real roots of this equation depends on 
the sign of the discriminant: 



which is positive in our case. That is why Eq. (33) has a real root. Real roots of Eq. (33) can 
be identified easily with the help of well-known Cardano's formula 22|. However, as it will 
become evident below, they are not necessary for a further specification of the expressions 
for the stationary points need for the study of the stability conditions. 

To address the question concerning the stability of the stationary points more precisely we 
linearize the set of Eqs. (29) in the vicinity of the stationary points Af^ = pq cosOq, A'^^ = 
Pq sinOo and Bq = 0, and obtain: 

5ii = -7(1 + 2rAf^Af^)6Ai - 7r(4°^' + 3Af ')M2, 

6A2 = 7r(Af ' + 3AS°^')Mi - 7(1 - 2rAf)Af )M2. (35) 
Alternatively, by taking into consideration the transformation (30): 

5A1 = RuSAi + R12SA2, 

6A2 = R21SA1 + R22SA2, (36) 

where 

^^f -7(1 + rpl sin 29o) -jrpl{3 - cos 29o) \ ^^^^ 
\ 7rp2(3 + cos 2^0) -7(1 - rpl sin 2^0) / ' 
As discussed in [l^, the type of the stability is determined by three characteristics of the 
matrix (R): 

T = Rii + R22, d = R11R22 ~ R12R21, T"^ — 4:d. (38) 

With the help of the matrix (/?), it is easy to check, that in our case the characteristics of 
the matrix are 

T = -27 < 0, d = 7^(1 + 8rVo) > 0, - Ad = 327VV0 > (39) 

and point towards the condition of a stable focus. Note, that the conditions (39) do not 
depend on the field parameter s = , and therefore hold for arbitrary values of the 
amplitudes {Ai, A2). 
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Thus, in the stationary resonance regime, when frequency of the external driving field is 
in the resonance with one of the modes of the two strongly coupled resonators, the stationary 
points are characterized by a stable focus. Therefore, we can argue that the dissipation leads 
to a stabilization of the secular instability regime. 



VI. ENERGY REDISTRIBUTION BETWEEN RESONATORS 

In this section we will address the problem of the energy redistribution between the 
resonators (Al/A^), which are pumped via the external fields. Let us suppose that the 
harmonic force acts only on the second resonator F2 = -F, Fi = 0. Then, with the help 
of the expression (10), for the ratio of the oscillation amplitudes we obtain the following 
relation: 



A, 



(40) 



We recall that the expressions (10) for the amplitudes Ai,2 with respect to the mode fre- 
quencies 1^1,2 have the same resonances embedded in the denominators. They naturally 
compensate each other and therefore do not appear in the ratio A1/A2. Nevertheless, as we 
see from Eq. (40), another resonance Wi ~ appears in the denominator of the expression 
Ai/A2. 

At first we neglect the influence of the damping and the nonlinearity terms, assuming 
that the frequency of the harmonic force is tuned with one of the partial frequencies of the 
resonators. Then for the case when F2 = F, fi ^ CJ2 and Fi = 0, we obtain 



A, 



A, 



D a 



(41) 



\ujI — 2 ' 

Hence, the relation between the amplitudes Ai and A2 is linear. From the second resonator 

2 

a fraction ^ of the energy is transferred to the first one. The damping and the nonlinear 
corrections can be considered again with the help of the substitution: 

^ , ^ B=^-^ (421 

CjI-^^ 2cI;iv/(/3A?- A)2 + 72' 4wi' ^ ' 

From now on we assume that 

71 ~ 72 = 7, xi = X2 = X and A = W2 - > 0. 
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FIG. 2: Energy redistribution curve between the coupled resonators and its asymptotics, plotted 
using Eq. (44) for the following values of parameter /3 = 0.51 • 10^*^ Hz/m, A = 0.17- 10*^ Hz, 7 = 
2.0 • 10^ Hz, f = 1.64 • 10^° Hz^. 

The resonant denominator is an important feature of the expression (42). When Ai changes 
the resonance condition holds in the expression (42). Performing the substitution (42) in 
the expression (40) and raising to the square, we obtain 

X f 



y {I3x - A)2 + 7' 



2 



(43) 



where x = A\, y = A^, f = So, instead of studying the dependence y = Y{x), for 
convenience one can convert the problem to one studying the following implicit function 

F(x, y) = x[{/3x - A)2 + 7^] -fy = 0. (44) 

By setting the derivative dy/dx equal to zero, one obtains an equation for the extrema of 
the function y = Y{x) 

fx = ^ /"'(3/3'^' - 4/3a:A + + f ) = 0. (45) 

It follows then that the points of extremum are 



2A/ / 3A2 + 72N 
14 



For simplicity, we consider the limiting case 7 ^ A. In this case we get two real roots from 
Eq. (46): 



'^^<0^ (48) 



It is easy to determine the signs of the second derivatives: 

fy =^>o ^ = 

dx"^ x=xi f dx'^ x=X2 f 

Therefore, the function y = Y{x) has a maximum at the point x = X2, and a minimum at 
X = Xi. The curve y = Y{x) is characterized by two asymptotes as well. The first one, 
for small values of x and y, is a linear function y = (A^//)x = (4/cr^)x. The second one, 
for large values of x and y, is a cubic function y = (y)x^. Using the results obtained in 
this section, one can plot the curve of the energy redistribution between the resonators, see 
Fig. 2. 

The anharmonicity of resonators' oscillations can significantly change the energy redis- 
tribution between the resonators. It turns out, that the energy pumped into the second res- 
onator via the external energy source F2 = F is transformed to the first resonator {Fi = 0) 
in a different way depending on the oscillations amplitude. For small amplitudes of the nor- 
mal modes, the energy transfer between the resonators is linear, Al = and the transfer 

rate is defined by the values of the connectivity a. With increasing oscillation amplitude the 

/ 2 N 1/3 

linear law is changed and turns nonlinear A\ = f | ■ -^A\ \ . Therefore, we can conclude 
that the anharmonicity of the oscillations degrades the energy transfer rate. 



VII. APPLICATION TO MASS MEASUREMENT SENSORS AND THE NON- 
LINEAR SHIFT OF THE MODE FREQUENCIES 



Nanomechanical resonators can be used as apprehensible sensors in many applications. 
For a review see 23|. A decisive advantage of the nanomechanical resonators are their 
resonance frequencies u ~ IGHz and quality factors Q ~ 10^ — 10^, which are significantly 
higher than those of electrical resonant circuits. That is the reason why nanomechanical 
resonators are sensitive transducers for the detection of molecular systems, in particular for 
biological molecules 16|. Resonant mass sensor devices operate by measuring the frequency 
shift which is proportional to the mass of the molecules of the material under investigation 
{2^. Details of the measurement protocol can be found in 



25|. Here we briefly refer to the 
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main facts. Assuming that the added mass 6M is smaller than the effective resonator mass 
M one can write a linearized expression 5M ^^5oj. The minimal measurable frequency 
shift 5uj naturally defines the sensitivity of the sensor. Due to thermal fluctuations 5uj > 0. 
For the sin gle, simple damped harmonic oscillator system the minimal measurable frequency 
shift reads 25 1 



iMu^A^ Q \ ■ ^ ^ 

Here A/ is the measurement bandwidth, M is the resonator mass, u is the frequency of the 
oscillation and T is the temperature. As follows from the analysis of the preceding sections, 
at low temperatures the nonlinear effects (that were not considered in 25|) can produce 
a frequency shift larger than the minimum measurable frequency shift associated with the 
thermal effects (see Eq. (18)) 61^2 > <^i,2- We propose to use the system of coupled nonlinear 
oscillators to act as an amplifier for the frequency shifts. We are convinced that in this way 



1/2 



far better mass measurements are possible for experiments described in 



25|. 



VIII. CONCLUSION 



We have developed a general analytical treatment of a system of two coupled driven non- 
linear nanomechanical resonators, which is valid for an arbitrary coupling strength (connec- 
tivity) between them. We derive general analytical expressions for the amplitude-frequency 
characteristics of the system with a special emphasis on the energy redistribution and the 
energy transport between the resonators. The obtained results are valid for arbitrary values 
of the connectivity. In the limit of a weak coupling one recovers the previously obtained re- 
sults jl^. In particular we have shown that the amplitude-frequency characteristic consists 
of two tilted peaks, the frequency separation between which is equal to twice the value of the 
resonators coupling constant 2D. If the frequency of the external force Q is swept the oscil- 
lation amplitude shows two hysteresis loops in the vicinity of the mode frequencies. These 
hysteresis loops contain unstable areas, in which a slight change of the driving frequency 
is accompanied by an instantaneous and a significant change of the oscillation amplitude. 
This is an interesting phenomenon, since it can be utilized to switch easily between the 
energy transport regimes of the resonators. We found that for small oscillation amplitudes 
the energy transfer between the resonators follows a linear law = jA\ and the transfer 
rate is entirely defined by the values of the connectivity a. With increasing the oscillation 
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amplitude the energy transfer law turns nonlinear Af = f | ■ ^^2) therefore the 

transport rate becomes slower. Switching off the energy transfer rate by tuning of the driv- 
ing field frequency is a simple protocol from an experimental point of view and therefore we 
expect it to be easily observable. 
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Appendix A — details from section IV 

By taking into account damping effects, the amplitudes of the forced oscillation of the 
nonlinear resonators for an arbitrary value of the connectivity a are: 
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Appendix B — details from section V 



The explicit form of the set of equations 
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